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LONG-TERM  GOALS 

This  proposal  consisted  of  six  individual,  but  interelated  projects.  A  common  component  of  all  tasks 
was  the  exploitation  of  satellite  imagery.  The  six  projects  are  listed  as  the  following  tasks: 

Task  1:  Air-Sea  Interactions  in  Typhoon  Conditions 
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Task  2:  Satellite  Observations  of  Surface  Signatures  and  Properties  of  Internal  Waves  in  Straits 

Task  3:  Satellite  Observations,  Surface  Signature  and  Properties  of  Nonlinear  Internal  Waves 

Task  4:  Satellite  Observations  of  Surface  Signatures  in  the  Philippine  Sea 

Task  5:  Satellite  Observations  of  Oceanographic  Processes  for  Quantifying,  Predicting,  and 
Exploiting  Uncertainty  (QPE) 

Task  6:  Exploitation  of  Commercial  Satellite  Radars 

The  long  term  goal  of  this  project  was  to  employ  satellite  SAR  imagery  to  the  quantitative  analysis  of 
typhoons,  internal  waves,  ocean  surface  features,  sea  state  prediction  and  coastal  dynamics.  To  achieve 
this,  a  variety  of  satellite  SAR  sensors  with  different  imaging  frequencies,  multi-polarizations  and 
different  beam  modes  were  used  to  determine  which  combination  of  the  above  capabilities  could 
provide  the  best  solution. 

OBJECTIVES 

For  simplicity  we  provide  only  a  primary  objective  of  each  task. 

a)  Task  1:  To  determine  the  characteristics  of  typhoons  with  satellite  SAR  imagery. 

b)  Task  2:  To  complement  the  acquisition  and  analysis  of  conventional  radar  intensity  images  and 
other  data  by  direct  high-resolution  imaging  of  surface  current  fields  from  space  by  along-track 
InSAR. 

c)  Task  3:  To  measure  surface  signatures  and  Bragg  modulations  during  the  passage  of  soliton 
packets. 

d)  Task  4:  To  improve  the  characterizations  of  the  Philippine  Archipelago  Straits  dynamics  we 
propose  to  acquire  satellite  data,  in  particular  synthetic  aperture  radar  (SAR)  imagery  of 
features  and  processes  present  in  these  straits. 

e)  Task  5:  To  determine  the  structure,  the  seasonal  and  annual  variability  of  the  cold  dome  over 
the  outer  continental  shelf  northeast  of  Taiwan  and  the  position  and  variability  of  the  Kuroshio 
Current,  in  particular  excursions  onto  the  continental  shelf. 

f)  Task  6:  To  explore  innovative  uses  for  commercial  radar  applications. 

APPROACH 

Satellite  based  Synthetic  Aperture  Radar  (SAR)  can  operate  day  and  night  and  in  all  weather 
condition.  This  is  of  importance  especially  in  the  Southeast  Asia  region  where  most  of  the  time  there  is 
a  persistent  cloud  cover  and  optical  sensors  such  as  MODIS  to  image  internal  waves  are  not  frequently 
useful.  However,  radars  penetrate  clouds  which  ultimately  are  useful  in  typhoons.  Furthermore,  radars 
are  truly  sensing  the  surface  and  therefore  observe  the  impact  of  the  atmosphere  on  the  ocean  surface. 
SAR  can  provide  large  spatial  coverage  up  to  500  km  wide  swaths  using  the  ScanS AR  mode  but  at 
lower  resolution  (~50  to  100  m)  for  the  directional  wave  spectra,  surface  wind  vectors,  and  surface 
roughness.  The  high-resolution  of  SAR  wind  retrievals  make  these  data  suitable  for  understanding 
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tropical  cyclone  morphology.  In  particular,  the  wind  structure  of  eye  wall  replacement  and  the 
formation  of  thunderstorms  and  tornados  in  hurricanes  can  be  imaged  from  all-weather  SARs. 

SAR  imaging  of  internal  waves  provides  the  necessary  geospatial  context  in  terms  of  internal  wave 
(IW)  characteristics  such  as  number  of  solitons,  wave  and  crest  length,  refraction  and  shoaling,  and 
origins.  These  images  were  also  used  to  guide  ship  and  glider  operations  in  areas  where  IWs  were 
present. 

Task  1:  From  July  to  October  2010  a  field  experiment  was  carried  out  east  of  Taiwan  in  the  primary 
region  with  high  probability  of  typhoon  tracks  (Figure  la).  As  part  of  the  experiment  we  deployed  two 
pairs  of  EASI  (Extreme  Air-Sea  Interaction)  and  ASIS  (Air-Sea  Interaction  Spar)  buoys  (Figure  lb)  in 
international  waters  about  200  km  apart.  The  buoys  were  equipped  with  extensive  meteorological  and 
oceanographic  instrumentation  to  provide  detailed  measurements  of  the  fluxes  and  properties  of  the 
air-sea  interface. 

We  also  utilized  two  remote  sensing  techniques  which  consisted  of  satellite-based  imagery  from 
RadarSat-2,  TerraSar-X,  Envisat  ASAR,  Cosmo-SkyMed  and  ERS-2;  and  a  ship-borne  marine  X-band 
system  (WaMoS)  on  the  R/V  Revelle.  Figure  2  shows  a  summary  of  the  quicklook  images  collected 
during  Typhoon  Megi.  Using  multiple  satellite  systems  allowed  for  a  good  coverage  of  the  storm  after 
entering  the  experimental  domain  and  across  the  Philippines  and  even  into  the  Taiwan  Straits.  The 
synthetic  aperture  radar  (SAR)  images  are  especially  valuable  in  determining  the  eye  location  and 
characteristics  of  storms.  The  radar  backscatter  data  is  used  to  derive  wind  speed  and  direction. 
Similarly  the  radar  data  can  be  used  to  estimate  wave  heights  and  wave  directions.  The  marine  radar 
systems  on  ships  can  also  derive  wind  speed  and  direction,  directional  wave  properties  including 
spectra  from  is  grazing  angle  observations.  The  range  of  measurements  from  ship  radar  is  more  limited 
and  extends  to  about  3  nm. 


Figure  1:  a)  Deployment  locations  of  the  EASI/ASIS  buoys  in  the  Western  Pacific  with  tracks  of 
three  of  the  typhoons  passing  by  during  the  2010  typhoon  season,  b)  Deployed  EASI/ASIS  buoys  in 

the  Western  Pacific. 
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Figure  2:  Multi  satellite  radar  collections  along  the  track  of  Typhoon  Megi. 

Task  2  and  Task  3:  Marine  radars  operate  by  transmitting  and  receiving  pulses  of  microwaves  that 
interact  with  the  rough  sea  surface.  The  radar  echo  from  the  sea  surface  is  primarily  dependent  on 
small-scale  ripple  waves  due  to  Bragg  scattering.  IWs  are  accompanied  by  a  spatially  and  temporally 
varying  surface  current  field.  Marine  radars  show  IWs  as  adjacent  bands  of  weak  and  strong 
backscatter  that  propagate  with  the  IW  packet’s  velocity.  As  for  SAR  images,  such  features  can  either 
result  from  hydrodynamic  interactions  between  the  short  surface  waves  and  the  horizontal  surface 
currents  caused  by  IWs,  or  from  surfactants  that  accumulate  in  IW  flow  convergence  zones  and 
dampen  the  short  surface  waves  there. 

A  novel  method  called  the  Radon  Transform  (RT)  was  applied  to  retrieve  IW  information  from  marine 
radar  data.  First,  the  images  were  divided  by  their  respective  radar  backscatter  ramps  to  produce  near- 
homogeneous  intensities.  The  ramp-corrected  images  were  then  used  to  retrieve  IW  packet  speed  and 
direction,  employing  a  combination  of  Radon  Transform  and  cross-correlation  techniques.  In  a  next 
step,  the  images  were  geo-referenced  and  corrected  for  the  IW  packet  motion.  By  averaging  a  4  minute 
sequence  of  such  preprocessed  radar  images,  we  can  now  significantly  enhance  the  IW  signatures.  The 
resulting  IW-enhanced  images  are  fed  into  our  IW  signature  retrieval  algorithm,  where  we  used 
localized  Radon  Transform  techniques  to  trace  the  signatures  of  each  individual  soliton,  yielding  IW 
soliton  maps.  From  such  maps,  one  can  easily  extract  parameters  such  as  wavelength,  velocity,  and 
curvature  for  each  IW  soliton. 

Figure  3  (left)  shows  a  transect  through  a  polar  X-band  radar  image  of  an  IW-enhanced  image  with  its 
IW  peaks  and  troughs  marked  by  black  and  red  dashed  lines,  respectively.  An  iterative  approach  that  is 
based  on  a  localized  Radon  Transform  technique  is  applied  to  track  each  soliton,  starting  from  the 
peaks  found  within  the  transect.  A  zoom  into  our  IW-enhanced  image  with  the  retrieved  IW  signatures 
as  well  as  the  position  of  the  transect  are  shown  in  Figure  3  (right). 
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Figure  3:  Transect  with  IW peaks  and  troughs  marked  by  dashed  lines  (left)  and  corresponding 
image  with  the  transect  position  marked  by  a  yellow  dashed  line  and  the  detected  IW  signatures 

marked  by  white  and  red  solid  lines  (right). 

Task  4:  The  Philippine  Archipelago  is  known  to  generate  strong  gap  jets  blowing  through  the  valley  of 
its  volcanic  mountains.  ENVISAT  synthetic  aperture  radar  (SAR)  images  during  the  2004-2009 
winter  monsoon  seasons  (November-March)  were  analyzed  using  the  wind  field  algorithm  derived  in 
Task  1.  The  wind  fields  revealed  the  presence  of  strong  jets  emanating  from  the  mountain  gaps  in  the 
lee  of  the  Philippine  Archipelago.  Jet  characteristics  exhibited  similar  features  of  classical  jets  in  water 
bodies.  The  Philippine  Island  gap  winds  were  compared  to  those  from  theoretical  literature  and 
experimental  studies.  For  all  jets,  the  initial  jet  width  dictated  the  distance  at  which  a  jet  became  self¬ 
similar,  wherein  a  larger  initial  jet  width  required  a  greater  distance  for  self-similarity.  Figure  4  shows 
examples  of  gap  jet  wind  fields  derived  from  the  SAR  backscatter  data.  The  jet  features,  long  streaks 
emanating  out  of  the  valleys  are  quite  obvious  in  contrast  to  the  low  wind  speeds  in  the  shadow  (lee)  of 
the  mountain  peaks. 
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Figure  4:  Examples  of  (a)  classic  (black  box)  and  decreasing  (white  box)  jet  cases  on  March  16, 
2009  1:58  UTC  and  (b)  increasing  jet  case  on  February  10,  2008  14:10  UTC.  The  background  is 
SAR-derived  wind  speed  (m/s),  with  vectors  denoting  SAR-derived  wind  direction.  Rectangular 
boxes  denote  the  jets  of  interest.  The  green  line  denotes  jet  half-width  (km),  the  perpendicular  black 
line  represents  initial  jet  width  (d,  km),  the  straight  black  line  is  the  jet  centerline,  and  the  blue  line 
represents  wind  speed  maxima  (m/s)  coordinates  associated  with  the  Principal  Component  Analysis 
technique  from  which  the  jet  centerline  is  fit.  Here  d=6  km  for  (a)  and  d=5  for  (b). 
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Task  5:  The  Ocean  northeast  of  Taiwan  has  a  large  continental  shelf,  complex  slope  bathymetry,  an 
upwelling  region,  a  variable  coastal  current,  a  strong  boundary  current  and  abundant  internal  waves. 
Intense  fishing  activity  inhibits  comprehensive  mooring  measurements.  These  features  make  it 
difficult  to  understand  the  physical  processes  of  this  region.  In  support  of  the  main  field  program  for 
the  QPE  DRI  in  August/September  2009,  15  long  strips  of  SAR  images  were  acquired  over  the 
continental  shelf  northeast  of  Taiwan  and  the  Kuroshio  Current  east  of  Taiwan.  The  collection  focused 
on  ascending  and  descending  passes  for  ERS-2  and  TerraSAR-X  data  using  different  beam  modes  that 
would  include  ScanSAR,  standard  and  fine  beam  modes  depending  on  the  location.  Similarly  we 
collected  MODIS  SST  images  and  other  EO  sensors  that  have  infra-red  channels  such  as  ASTER 
(Advanced  Spacebome  Thermal  Emission  and  Reflection  Radiometer),  a  high  resolution  imaging 
instrument  flying  on  Terra. 

Figure  5  shows  a  TerraSAR-X  image  over  the  northeast  Taiwan  shelf  with  abundant  internal  waves 
and  soliton  interaction  patterns.  Non-linear  internal  waves  are  observed  in  the  satellite  imagery  and 
are  used  to  deduce  propagation  patterns.  The  internal  waves  in  this  region  are  characterized  by 
elaborate  wave  patterns  created  by  multiple  internal  wave  propagation  directions  and  intersection 
angles.  Although  the  images  acquired  show  complicated  internal  wave  patterns,  three  preferred 
pathways  were  found  to  exist.  The  first  is  offshore  to  the  southeast  toward  deep  water.  The  second  is 
northwest  toward  the  coast  of  China  similar  to  internal  wave  patterns  observed  in  other  shelf  regions. 
The  third  is  orthogonal  to  the  first  two,  propagating  to  the  northeast  along  the  shelf  break.  This  third 
pathway  produces  the  complex  patterns  frequently  observed  in  satellite  imagery  of  the  region. 

25*30' 


25*15’ 


25*00' 

122*06’  122*12'  122*18’ 

Figure  5:  A  TerraSAR-X  image  collected  on  9  August  2009  showing  interacting  internal  wave 
patterns  and  soliton  groups  propagating  in  multiple  directions. 
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The  majority  of  internal  waves  along  the  shelfbreak  are  classic  depression  waves  generated  by  the 
internal  tide  along  the  shelfbreak.  The  remaining  internal  wave  signatures  are  possibly  generated  by 
the  encroachment  of  the  Kuroshio  Current  on  the  shelf  or  the  persistent  cold  dome.  Additional  internal 
waves  may  be  trans-basin  waves  advecting  onto  the  shelf  from  the  Ilan  Ridge  or  may  be  remnants  of 
interacting  wave  packets  (Figure  6).  From  the  analysis  of  the  imagery  we  found  220  internal  wave  (IW) 
packets  with  more  than  2,407  waves  (cf.  Figure  6).  There  are  clearly  some  hotspots  of  persistent 
internal  waves  located  on  the  shelf  and  over  the  entrance  to  the  canyons. 


Figure  6:  Locations  of  internal  wave  fronts  identified  from  SAR  imagery  acquired  between  2007 

and  2010  from  ERS-2  and  TerraSAR-X. 


The  histogram  of  propagation  directions  (Figure  7)  shows  the  primary  pathways  are  generally  oriented 
perpendicular  to  the  bathymetry.  The  majority  of  the  internal  wave  packets  are  traveling  southeast  (90° 
- 150°)  toward  the  open  ocean.  The  secondary  propagation  direction  is  to  the  northeast  (270°  -  330°) 
toward  the  coast.  This  propagation  direction  is  frequently  seen  along  the  shelf  in  the  South  China  Sea 
and  the  New  York/Mid  Atlantic  Bight.  A  tertiary  pathway  exists  towards  the  northeast  (30°  -  90°) 
along  the  shelfbreak.  Note  that  the  range  of  propagation  directions  along  shore  and  offshore  results  in 
a  peak  in  the  histogram  at  the  arbitrary  delineation  of  90°  toward  the  east. 
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Figure  7:  The  propagation  direction  of  internal  waves  from  historical  data.  The  northwest  or 
onshore  direction  and  southeast  or  offshore  directions  are  dominant.  However,  internal  waves  also 
travel  to  the  southwest  toward  the  east  coast  of  Taiwan.  Internal  waves  also  propagate  toward  the 

northeast  along  the  shelfbreak. 


RESULTS 

This  projects  generated  13  peer-reviewed  publications  and  4  more  in  review  as  well  as  5  conference 
proceedings.  In  addition  the  project  supported  four  graduate  students  of  which  three  received  their 
PhDs  and  one  more  will  finish  in  May  2015.  The  project  also  supported  two  post-doctoral  scholars. 

Below  is  a  list  of  doctoral  thesis  and  journal  article  publications  loosely  associated  with  each  task.  We 
also  include  the  abstract  for  each  publication. 

DOCTORAL  THESIS 

Bjorn  Lund,  October  2012,  Principal  Advisor,  (RSMAS).  Thesis:  Currents  from  X-Band  Marine  Radar 
Measurements  of  Internal  Waves 

Henry  Potter,  March  2014,  Committee  Member,  (RSMAS).  Thesis:  A  Study  of  Turbulent  Processes  at 
the  Air-Sea  Interface  in  High  Wind  Speeds 

Clarence  O.  Collins,  III,  August  2014,  Principal  Advisor,  (RSMAS).  Thesis:  Typhoon  Generated 
Surface  Gravity  Waves  Measured  by  Nomad-Type  Buoys 

Astrid  Werkmeister,  expected  May  2015,  Principal  Advisor,  (RSMAS).  Thesis:  Vertical  Structure  and 
Wind-  Wave  Interaction  in  the  Boundary  Layer  of  Tropical  Cyclones  Derived  from  SAR  Imagery 
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typhoons  on  the  ocean  in  the  pacific:  ITOP.  Bull.  Americ.  Meteo.  Soc.,  95(9),  1405-1418.  doi: 
http://dx.doi.Org/10.1175/BAMS-D-12-00104.l 
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2014:  EASI:  An  air-sea  interaction  buoy  for  high  winds.  J.  Atmos.  Oceanic  Technol.,  31,  1397- 
1409.  DOI:  10.1 175/JTECH-D-13-00201.1 

Romeiser,  R.,  H.C.  Graber,  M.J.  Caruso,  R.E.  Jensen,  D.T.  Walker,  A.T.  Cox,  2014:  A  New 

Approach  to  Ocean  Wave  Parameter  Estimates  From  C-Band  ScanSAR  Images.  IEEE  Trans. 
Geosci.  Remote  Sens.,  53(3),  1320  -  1345. 
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